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Abstract 

We report the first observation of the charmless hyperonic B decay, — > pA-Tr", using a 78 fb^^ 
data sample recorded on the T(4S) resonance with the Belle detector at KEKB. The measured 
branching fraction is B{B^ pAtt") = (3.97l^:^[J ± 0.56) x 10"^ Searches for B° pAK' 
and pS°7r~ yield no significant signals and we set 90% confidence-level upper limits of B{B'^ 
pAK-) < 8.2 X lO-'^ and B{B^ pfP-K') < 3.8 x 10"^ 

PACS numbers: PACS numbers: 13.20.H 
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The Belle collaboration recently reported the observation of ppK^ which is 

the first known example of B meson decay to charmless final states containing baryons. The 
three-body decay rate is larger than the rate for two-body decays (such as 5 — >• 0), 
and the observed Mpp spectrum peaks near threshold, in agreement with theoretical sugges- 
tions 0, 1^ . In this Letter we report the first observation of the related three-body decay 

— >■ pK'K~ , and a search for B^ pAK~ and pS°7r~ modes. The rate for 5° pAvr" is 
comparable to B~^ ppK~^, and the observed M^j^^ spectrum again peaks toward threshold. 

In the Standard Model, these decays proceed via h ^ u tree and h — >• s{d) penguin 
diagrams. They may be used to search for direct CP violation and test our theoretical 
understanding of rare decay processes involving baryons 0, B IE IE 0]- Modes involving 
hyperons, in particular, can probe the s quark chirality in B decay with sufficient 
statistics, they could provide a tool for probing T violation via the self-analyzed A 
polarization information. 

We use a 78 fb~^ data sample, consisting of 85.0 ± 0.5 million BB pairs, collected by 
the Belle detector at the KEKB asymmetric energy e^e~ (3.5 on 8 GeV) collider. The 
Belle detector is a large-solid-angle magnetic spectrometer that consists of a three-layer 
silicon vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel 
threshold Cerenkov counters (ACC), a barrel-like arrangement of time-of- flight scintillation 
counters (TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals located 
inside a super-conducting solenoid coil that provides a 1.5 T magnetic field. An iron fiux- 
return located outside of the coil is instrumented to detect mesons and to identify muons. 
The detector is described in detail elsewhere 0. 

Since the e^e~ center-of-mass energy is set to match the T(4S) resonance, which decays 
into a BB pair, one can use the following two kinematic variables to identify the recon- 
structed B meson candidates [lE| : the beam-energy constrained mass, Mbc = \J E"^^^^ — 
and the energy difference, = Eb — -Ebeam, where -Ebeam is the beam energy, and pb and 
Eb are the momentum and energy of the reconstructed B meson in the T(4S) rest frame. 
The candidate region is defined as 5.2 GeV/c^ < Mbc < 5.29 GeV/c^ and |AE| < 0.2 GeV 
in this analysis. 

The event selection criteria are based on the information obtained from the tracking 
system (SVD-I-CDC) and the hadron identification system (CDC-I-ACC+TOF), and are 
optimized using Monte Carlo (MC) simulated event samples. 

All primary charged tracks are required to satisfy track quality criteria based on the track 
impact parameters relative to the interaction point (IP). The deviations from the IP position 
are required to be within ±0.3 cm in the transverse {x-y) plane, and within ±3 cm in the z 
direction, where the z axis is defined by the positron beam line. Primary proton candidates 
are selected based on p/K/ir likelihood functions obtained from the hadron identification 
system. We require Lp/^Lp + Lx) > 0.3 and Lp/ {Lp + L^,) > 0.6, where Lp/K/n stands for the 
proton/kaon/pion likelihood. For kaons (pions), we require the kaon (pion) K-n likelihood 
ratio to be greater than 0.6. A candidates are reconstructed via the pvr" decay channel using 
the method described in Ref. 0. E° candidates are reconstructed via the A7 decay channel, 
where we use a 35 MeV/c^ mass window around the nominal mass and require the 7 
energy to be greater than 100 MeV. 

The dominant background for the rare decay modes reported here is from the continuum 
e~^e~ — s> qq process. The background from B decays is much smaller. This is confirmed 
with an off-resonance data set (8.8 fb~^) accumulated at an energy that is 60 MeV below 
the T(4S), and an MC sample of 120 million continuum events. In the T(4S) rest frame. 
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continuum events are jet-like while BB events are spherical. We follow the scheme defined 
in Ref. 3| and combine seven shape variables to form a Fisher discriminant fisl] that is 
used to optimize continuum background suppression. The variables chosen have almost 
no correlation with Mbc and AE. Probability density functions (PDF's) for the Fisher 
discriminant and the cosine of the angle between the B flight direction and the e~ beam 
direction in the T(4S) rest frame are combined to form the signal (background) likelihood 
'^s(BG)- We require the likelihood ratio CTZ = C^/i'^s + -^bg) to be greater than 0.8, which 
suppresses about 94% of the background while retaining 66% of the signal. The signal and 
background PDF's are obtained from MC simulation studies. 
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FIG. 1: The (a) AE and (b) Mbc distributions for pK-n~ candidates. The solid, dotted and 

dashed lines represent the combined fit result, fitted background and fitted signal, respectively. 



Figure 1(a) shows the AE distribution for selected pAvr candidates that have Mbc > 
5.27 GeV/c^; Fig.I^b) shows the Mbc distribution for events with \ AE\ < 0.03 GeV. With 
the current statistics, no intermediate resonances are evident in the Dalitz plot for this 
channel. We use a binned likelihood flt to estimate the signal yield. A Gaussian is used 
to parameterize the signal in Mbc while a double Gaussian is used for AE. The Gaussian 
parameters are determined from MC simulation. Background shapes are studied using side- 
band events (0.1 GeV < \AE\ < 0.2 GeV for the Mbc study and 5.20 GeV/c^ < Mbc < 5.26 
GeV/c^ for AE ) and are checked with the continuum MC sample. We use the ARGUS 
function [l3| to model the Mbc background and a linear function for the AE background. 
The fit results are shown as curves in Fig. ^ The fit to the AE distribution yields 39.2 
candidates with a significance of 5.8 standard deviations. The fit to the Mbc distribution 
yields 33.7^7^ candidates with a significance of 5.7 standard deviations. The smaller yield in 
the Mbc fit is consistent with the AE fit result after taking into account the efficiency of the 
|Aii^| < 0.03 GeV selection. The signal yields and the branching fractions are determined 
from fits to the AE distribution rather than to Mbc in order to minimize possible bias from 
BB background, which tends to peak in Mbc but not in AE. 
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Since the decay is not uniform in phase space, we fit the Ai? signal yield in bins of M^^, 
and correct for the MC-determined detection efficiency for each bin. This reduces the model 
dependence of the branching fraction determination. The signal yield as a function of pK 
mass is shown in Fig. |21 The distribution from a three-body phase space MC, normalized 
to the area of the signal, is superimposed. The observed mass distribution peaks at low pK 
mass, similar to that observed for — >■ ppK~^ decays The results of the fits, along with 
the efficiencies and partial branching fractions for each M^;^ bin, are given in Table HI We 
sum the partial branching fractions in Table |l] to obtain 



B{B^ pAn-) = (3.97^^:°° (stat) ± 0.56 (syst)) x 10" 
where the systematic uncertainty is described below. 



TABLE I: The event yield, efficiency, and branching fraction {B) for each bin 



M X ( GeV/c2) signal yield efficiency (%) B (10" 
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FIG. 2: The fitted yield divided by the bin size for — > pAvr 
distribution is from a phase-space MC simulation with area normalized to signal yield 



as a function of M^^. The shaded 
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The systematic uncertainty in particle selection is studied mainly using high statistics 
control samples. Proton identification is studied with a A — > pvr" sample. Kaon/pion 
identification is studied with a D*^ D^vr"*", K~tt^ sample. Tracking efficiency is 

studied with ^ 77 and 77 tt+tt^tt^ samples. Based on these studies, we assign a 2% error 
for each track, 3% for each proton identification requirement, and 2% for each kaon/pion 
identification requirement. 

We study the CTZ continuum suppression by varying the CIZ cut value from to 0.9 
to check the systematic trend. The systematic error is found to be 4%. The additional 
uncertainty of off-IP tracks for A reconstruction is estimated to be 6%, which is determined 
from the difference of the proper decay time distributions for data and MC simulation. 

The systematic uncertainty in the fit yield is studied by varying the parameters of the 
signal and background PDF's. We assign an error of 3% for this. The MC statistical 
uncertainty and modeling with eight M^j^ bins contributes a 4% error in the branching 
fraction determination. The error on the number of total BB pairs is determined to be 1%. 
The error of the branching fraction for A pTT~ is 0.8% [llj. 
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FIG. 3: The (a) AE and (b) M^c distributions for 
represents the fit results. 



pAK candidates. The solid line 



The tracking systematic error is estimated to be 8% by summing the correlated errors of 
2% per charged track. The particle identification error is estimated to be 8% by summing the 
correlated errors of 3% per proton identification and 2% for the primary pion identification. 
Then we combine them in quadrature along with other uncorrelated errors to determine a 
total systematic error of 14%. 

We also search for the decay modes S° — > pAK~ and — > pE^Tr". The Mbc and 
AE distributions with fit projections are shown in Figs. El and ^ With the signal region 
extended to \AE\ < 0.04 GeV and Mbc > 5.27 GeV/c^, no significant signal is found. We 
use the fit results to estimate the expected background, and compare this with the observed 
number of events in the signal region in order to set the upper limit on the yield at the 90% 
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FIG. 4: The (a) AE and (b) Mbc distributions for — > pT,^7r~ candidates. Note that events from 
pAiT^ decay can feed into the high AE region (~ 0.1 GeV). The corresponding distribution 
is obtained from the pAvr" MC and included in the fit. 



confidence level jl^, [l^ E3|- This procedure takes systematic uncertainties into account. 
The estimated backgrounds are 28.9 ± 2.6 and 50.5 ± 4.0, the numbers of observed events 
are 26 and 56, the systematic uncertainties are 14% and 28%, and the upper limit yields are 
8.3 and 22.4 for pAK~ and pS°7r~, respectively. We estimate the efficiencies from a phase 
space MC sample. The 90% confidence-level upper limits for the branching fractions are 
^ pAK-) < 8.2 X 10"^ and ^ pS^vr-) < 3.8 x 10"^ 

Following our observation of the ppK^ mode, some authors 0| predicted a 

much smaller branching fraction (< 10^^) for the B^ pA7i~ mode, but a relatively sizable 
B^ —>■ pTPti~ . Although the predicted rates are not borne out by our present findings, the 
threshold peaking behavior shown in Fig. |21was anticipated 0, 0, 0- 

In summary, we have performed a search for the rare baryonic decays B^ —>■ pAn^ , pAK~, 
and pE°7r~ with 85.0 ± 0.5 million BB events. A clear signal is seen in the pAvr" mode, and 
we measure a branching fraction of B{B^ pAn-) = (3.97lo:8o ± 0.56) x 10"^ The other 
two modes are not seen, and we set 90% confidence-level upper limits of B{B^ — > pAK~) < 
8.2 X 10"^ and B{B^ pE^ir-) < 3.8 x 10"^ 
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